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Preface 


The Soil and Water Management & Crop Nutrition (SWMCN) Subprogramme 
of the Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture 
supports IAEA Member States and FAO Member Countries in the use of fallout 
radionuclides (FRNs) techniques for assessing soil redistribution patterns and 
magnitudes to reduce the impact of soil erosion on agro-ecosystems. 

From all available FRNs (e.g. '°’Cs, 7!°Pb.x, 73°*?*°Pu), Beryllium-7 (Be) is the 
only soil tracer that permits the gathering of short-term soil redistribution infor- 
mation. This natural cosmogenic isotope was used in this context for the first time at 
the end of the 1990s. Developments over the past decades now permit investigation 
of erosion processes not only during short extreme weather events but also over 
extended periods of up to several months. 

The SWMCN Subprogramme has actively sought to refine the use of ‘Be in 
close collaboration with the School of Geography, Earth and Environmental 
Sciences at the University of Plymouth, UK, and with the help of scientists from 
developed and developing countries. This book reflects the latest developments 
obtained and the state of the art related to the use of this natural tracer of soil and 
sediment. 

This handbook presents the foundation of the Be method and provides guide- 
lines for the practical applications of basic and advanced approaches of this tech- 
nique. It assists scientists, technicians and students to assess the impact of soil 
erosion events and further support the development of land management and 
climate-smart agriculture policies for future sustainable food security. 


Vienna, Austria Lionel Mabit 
Plymouth, UK William Blake 
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Chapter 1 A) 
The Use of Be-7 as a Soil and Sediment E 
Tracer 


A. Taylor, W. H. Blake, A. R. Iurian, G. E. Millward and L. Mabit 


1.1 Origin and Agro-Environmental Behaviour of "Be 


1.1.1 Atmospheric Production of” Be 


Beryllium-7 (Be) (T12 = 53.3 days) is a cosmogenic fallout radionuclide (FRN) 
produced in the upper atmosphere by cosmic ray spallation of nitrogen and oxygen. 
As reported by Kaste et al. (2002), rates of production are dependent upon solar 
activity with greater production in the stratosphere than the troposphere and generally 
increased production occurring at higher latitudes owing to cosmic ray deflection 
towards the poles. 

Following its formation, "Be becomes associated with aerosols and its flux is 
then controlled by complex atmospheric transport processes, which display seasonal 
variations largely linked to atmospheric mixing and precipitation patterns. During 
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the spring months at mid latitudes, increased stratosphere-troposphere exchange can 
occur during folding of the tropopause leading to higher concentrations of ‘Be in the 
troposphere (Kaste et al. 2002). Once in the troposphere, ‘Be is subjected to vertical 
mixing through convective circulation, which is likely to increase during warmer 
months and serves to transport ‘Be enriched air to the lower troposphere, increasing 
its availability for precipitation scavenging. Given that precipitation scavenging is 
the main mechanism for the removal of ‘Be from the atmosphere, seasonal and lati- 
tudinal climatic conditions exert a strong influence upon atmospheric concentrations 
(Kusmierczyk-Michulec et al. 2015). 


1.1.2 7Be Fallout 


Wet deposition governs ’Be delivery to the Earth’s surface while dry deposition only 
accounts for approximately 10% of the overall fallout (Kaste et al. 2002). Activity 
concentration of 7Be in rainwater (Bq L7!) is conditional on its availability for 
scavenging and is influenced by a number of factors including rainfall rates and 
volumes (Ioannidou and Papastefanou 2006). 

The term ‘washout’ refers to situations when rainfall scavenging exceeds the 
availability of 7Be in the atmosphere leading to a reduction in Be rainwater activity 
concentration across an event and following high magnitude or prolonged rainfall 
events. In contrast, greater activity concentrations during low magnitude events can 
be related to efficient scavenging of "Be by fine droplets where there is an abun- 
dance of ’Be-bearing aerosols (Ioannidou and Papastefanou 2006). The processes of 
atmospheric circulation and depletion not only affect the activity concentrations of 
7Be on a temporal scale, but are also likely to increase its spatial variability (Taylor 
et al. 2016). "Be deposition (Bq m~’) is significantly correlated to the amount of 
rainfall received on an event basis and, therefore, seasonal depositional fluxes reflect 
monthly rainfall volumes (Doering and Akber 2008). 


1.1.3 7Be Sorption Behaviour in Soils and Sediments 


Upon deposition and infiltration into the soil surface, 7Be is known to be rapidly 
adsorbed to fine soil particles (Fig. 1.1), which is supported by laboratory batch 
studies, high partition coefficients (Kg > 10°) in aquatic systems and shallow depth 
distributions displayed in a range of soil types (Taylor et al. 2013). 7Be shows prefer- 
ential adsorption to fine sediment fractions (Taylor et al. 2014), which are likely to be 
readily mobilised in hillslope systems. In solutions of pH 5-6 and in the absence of 
humic acid, Be is likely to be present as Be?* and the hydrolysed species, BeOH*; 
both of these forms are highly reactive, resulting in rapid sorption rates (Kaste et al. 
2002). The rapid adsorption together with a short half-life (the latter precludes the 
effects of diffusion and bioturbation over time) leads to distributions typically extend- 
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Fig. 1.1 Schematic diagram of ’Be production and fallout. ’Be is produced by cosmic ray spallation 
in the upper atmosphere and becomes associated with aerosols. ’Be-bearing aerosols are scavenged 
by precipitation and, consequently, wet deposition is the dominant pathway to the Earth’s surface. 
7Be is rapidly adsorbed to soil particles and remains cohesive under oxic field conditions 


ing to around 20 mm below the soil surface and an exponential decrease with depth 
(Fig. 1.2). Where soils remain aerated and free from waterlogged conditions, as is 
the case in most agricultural field sites, Be is likely to remain adsorbed to soil par- 
ticles (Taylor et al. 2012). Both characteristics (i.e. rapid adsorption and cohesive 
behaviour) are key prerequisites for ‘Be application as a soil redistribution tracer. 


1.2 The Use of ’Be as a Soil Redistribution Tracer 


The application of FRN tracers (e.g. "Be, '°7Cs, *!°Pb,,, Pu) provides distinct and 
additional advantages over traditional soil monitoring techniques by enabling ret- 
rospective estimates of soil redistribution from relatively few site visits. Through 
carefully planned sampling programmes, high spatial resolution estimates of soil 
redistribution can be acquired, which are unlikely to be achieved using conven- 
tional methods (Walling et al. 1999). In this regard, the use of nuclear and iso- 
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Fig. 1.2 Example of 7Be 7Be activity (Bq kg") 
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topic techniques enables rigorous assessment of soil conservation measures and, if 
applied over extended time periods, can help to determine the sustainability of agri- 
cultural systems. 

FRN tracing techniques have been commonly applied to assess soil redistribution 
at the hillslope scale using !7Cs (T 1). = 30.2 years) and 7!°Pb (T 1. = 22.3 years). The 
half-life of these radionuclides enables redistribution estimates relating to medium- 
term (i.e. decadal) timescales, providing important information with regard to soil 
loss or gain in an historical context. However, the application of these decadal-scale 
radionuclides does not allow for determination of soil redistribution as a result of 
recent changes in land use or recent rainfall events. This is an important consider- 
ation given the need for assessment of soil conservation measures with increasing 
recognition of soil as a finite global resource. The short-lived 7Be is, therefore, com- 
plementary to the other FRN tracers in enabling estimates of soil redistribution across 
much shorter timescales such as rainfall events or wet seasons. 

The premise for using ‘Be as a soil redistribution tracer is that the tracer will be 
rapidly adsorbed to soil particles upon fallout, thus establishing an inventory ("Be 
activity per unit area) (Bq m~?) at a sampling site. A key assumption of the technique 
is that ’Be will remain adsorbed to soil particles across the timescale of the study and, 
therefore, any changes in inventory (i.e. depletion or gain) at a sampling point can 
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only occur as a result of fallout, radioactive decay and/or soil redistribution processes. 
The Be inventory at a site, which experiences soil redistribution, can be compared 
directly to the "Be inventory at a stable location, referred to as a reference site, where 
there is no loss or gain in inventory as a result of soil redistribution. The inventory at 
a reference site is established only through 7Be atmospheric fallout and, given that 
this is predominantly delivered by precipitation, reference sites are typically level 
areas which retain the fallout without loss except for radioactive decay. Alongside the 
important tracer assumptions of rapid adsorption upon fallout and cohesive behaviour 
during the period of study, it is also crucial that the reference site receives the same 
fallout deposition as the eroding area of interest, which requires both sites to be in 
close proximity. In addition to this, it is also assumed that the fallout received at 
the study site will be uniform so any differences in inventory across the study area 
can be attributed to soil redistribution. If these assumptions hold, then any reduction 
in inventory at a sampling location with respect to the reference inventory can be 
attributed to soil loss and a gain in inventory would suggest soil deposition, given that 
the redistribution of "Be (as a cohesive tracer), will only occur in line with movement 
of the sediment to which it is adsorbed. 

To quantify the mass of soil eroded or deposited (kg m~), it is essential to deter- 
mine the depth distribution of ’Be in the soil profile and, in particular, the relaxation 
mass depth (referred to as ho) a value, which describes the shape of the ’Be distri- 
bution. This is an essential component of the commonly applied soil redistribution 
model described below. On the basis of this, any change in ‘Be inventory at a sam- 
pling location can be fitted to the mass of soil redistribution required to implement 
that change. An example is given in Fig. 1.3, which shows a loss of 100 Bq m~? at 
the eroded location (i.e. inventory of 300 Bq m~?) with respect to the stable reference 
inventory (i.e. 400 Bq m~”). By establishing the exponential shape of the ’Be profile, 
the conversion model is able to estimate the mass depth of soil required to reduce 
the inventory by 100 Bq m7. 

The specific conversion model used to derive soil redistribution magnitudes (Pro- 
file Distribution Model [PDM]) is detailed in Blake et al. (1999) and Walling et al. 
(2009) and the basis of this is outlined in Appendix 1.1. In addition, a worked example 
of this model is provided in Chap. 4. 


1.3 Examples of the Application of ’Be as a Soil 
Redistribution Tracer 


Blake et al. (1999) and Walling et al. (1999) carried out the first comprehensive 
investigation into the use of ’Be to estimate short-term soil redistribution by studying 
the impact of a heavy rainfall event at an agricultural site in southwest England. The 
7Be depth distributions from this study site are shown in Fig. 1.4. The key requirement 
of uniform "Be fallout across the study location was satisfied owing to the soil Be 
inventory being reset (i.e. reduced to below detectable limits) by ploughing prior to 
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Fig. 1.3 Conceptual diagram of the premise for using Be as a soil erosion tracer. The figure shows 
the 7Be inventory (400 Bq m~?) and depth distribution at a stable reference location, which has 
been established by 7Be atmospheric fallout. The eroded and depositional zones (representing a 
neighbouring hillslope area for example) have received the same fallout as the reference location, 
thus, an inventory has been established but the zones have experienced depletion or enhancement 
of the 7Be inventory as a result of soil erosion and soil deposition respectively. This simplified 
conceptual diagram can be compared to examples of measured 7Be inventories in Fig. 1.4 


the commencement of the investigation. There was also a period of low intensity 
rainfall prior to the erosive event, which established a measurable Be inventory 
without causing significant soil redistribution. A reference baseline inventory was 
estimated from neighbouring undisturbed pasture locations against which the slope 
inventories could be compared. The depth profile was determined in a non-eroded 
part of the study field. The PDM was used to determine soil redistribution rates and 
results indicated high levels of soil export from the site as a consequence of the 
intense rainfall event and compacted soil conditions, highlighting the pressing need 
for soil conservation measures to maintain soil fertility and to prevent sediment from 
entering catchment watercourses. 

As highlighted by Mabit et al. (2014), an important advantage in the use of Be 
as a soil erosion tracer is the ability to assess the effectiveness of recent soil conser- 
vation measures or the impacts of land use change. Schuller et al. (2006) applied the 
approach proposed by Blake et al. (1999) to document the effects of recent forest 
clearance upon slope soils in a timber harvesting region of Chile. Through the ’Be 
approach, the use of mitigation strategies could be evaluated and, in this case, the 
need for further improvements in soil conservation strategies was identified. Sepul- 
veda et al. (2008) aimed to assess the effects of land management practice (in this 
case stubble burning) by using a combined FRN approach ("Be in conjunction with 
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Fig. 1.4 7Be depth distributions at a stable reference location (a), eroded hillslope locations (b and 
c) and a depositional hillslope location (d) [based on Walling et al. (1999)]. Note the depleted 7Be 
inventories at the eroding sites and the increased inventory at the depositional site with respect to 
the reference location 
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137Cs), Here, the combined use of FRNs was able to compare impacts associated 
with different land management scenarios across a range of timescales, providing 
key information with regard to the sustainability of the agricultural system. Blake 
et al. (2009) also applied a combined approach and used ’Be in conjunction with !37Cs 
and 7!°Pb to determine post wildfire sediment budgets at hillslope sites in southeast 
Australia. Data identified major sediment source areas and focussed attention on the 
implications of wildfire for water quality management. 

Blake et al. (1999) and Sepulveda et al. (2008) identified the occurrence of extreme 
erosion events by comparing short term erosion rates with medium-term estimates 
derived from !?’Cs. Walling et al. (1999) showed that medium-term erosion rates at 
the study site were approximately 5 times lower than the short-term rates, suggesting 
that the studied rainfall event may not have been representative of seasonal patterns 
in rainfall. Owing to this, adjustments to the model were suggested by Walling et al. 
(2009) to enable the PDM to be applied across a greater timescale, encompassing a 
number of erosive rainfall events and improving representativeness. Thus, a method 
was developed for extending the use of the PDM to periods covering a few months 
or wet seasons. Within this extended timescale approach, the fundamental aspects of 
the model remain similar to the conventional approach whereby ‘Be inventories at a 
sample point are compared to a reference baseline and the depth distribution of ’Be 
is measured. However, for the approach to be applied across multiple erosive rainfall 
events it is crucial to consider changes in the Be inventory across the study period 
to avoid underestimation of erosion rates. Additional components of the model are, 
therefore, included whereby loss of inventory at a sampling point is accounted for by 
estimating the relative erosivity of each rainfall event and ongoing radioactive decay. 
Similarly, assessment of inventory gain accounts for additional ‘Be fallout during the 
study period alongside soil redistribution. The application of this modified approach 
was piloted by Schuller et al. (2010), with soil conservation measures shown to be 
less effective across an extended time period, which encompassed a range of rainfall 
events, demonstrating the value of considering soil erosion in a long-term context. 

Recent plot studies have validated the method on a site-specific basis by demon- 
strating comparability between the ’Be-derived soil redistribution budget and direct 
measurements (e.g. Porto and Walling 2014; Shi et al. 2011). 


1.4 The Requirement for a Standardised Approach 


FRN tracers are likely to be used in a wide range of environments to assess soil 
redistribution and, as such, numerous methodological approaches for establishing the 
key model parameters have been applied. Although it is necessary to tailor sampling 
approaches to local environmental conditions, it is also important that the approach 
is standardised, as far as is practicable, to ensure data quality and aid comparability 
(Mabit et al. 2008). 

7Be has the potential to be a valuable decision support tool for catchment managers 
although to achieve robust estimates of soil redistribution using the PDM it is crucial 
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to accurately estimate its key parameters namely:(1) the reference ’Be inventory, 
and (2) the "Be depth distribution (the latter to ensure an accurate value of ho). It is 
also important to consider whether the selective transport of particle size fractions 
is occurring during rainfall events given preferential association of ‘Be to fine soil 
fractions (<63 um). Failure to account for preferential adsorption to fine particles 
can lead to overestimation of soil erosion rates since a large proportion of the Be 
inventory is associated with the fine sediment fractions, which are likely to be readily 
mobilised during rainfall events. Hence a relatively large deficit in "Be inventory 
could actually equate to a low mass of soil loss. 

The requirement for standardising procedures is of particular importance for deter- 
mining ’Be depth distributions in light of recent studies, which demonstrate how the 
choice of sampling method can influence results (Baumgart et al. 2016; Ryken et al. 
2016). Determining ’Be depth distributions in soils is a practical challenge and one 
that requires careful consideration of the sampling site, both in terms of stability 
and uniformity (Ryken et al. 2016) as well as a precise sampling method. The latter 
requires a technique that allows fine depth increments (2 mm depths) to be collected 
throughout the Be profile enabling enough detail to be able to accurately determine 
ho. Mabit et al. (2014) designed and tested a reproducible tool and a modus operandi, 
which permits collection of fine soil increments in a range of soil types and provides 
a basis for standardising the sampling method for ho determination. 

In practice, developing a robust sampling strategy to establish model parameters 
and assess soil redistribution is likely to require a balance between the number of 
samples required to ensure data quality and the capacity for sample analysis. As with 
any field operations, preliminary knowledge of the agro-environmental conditions 
and background data (e.g. land use history; soil hydrological properties) at the sam- 
pling site is beneficial and will help to streamline a pragmatic sampling programme. 

Against this background, Chap. 2 provides a guideline for undertaking a soil redis- 
tribution study using ‘Be as a tracer by outlining a standard approach for establishing 
key conversion model parameters. 


Glossary 

FRN Fallout Radionuclide. 

Becquerel The SI unit of radioactivity. 

Activity Concentration Radioactivity per unit of volume or mass e.g. Bq 
L-!; Bq kg™!. 

Inventory Radioactivity per unit area (Bq m~?) also termed 
areal activity. 

Relaxation mass depth (ho) The soil mass depth (kg m~?) at above which 
63.2% of the beryllium-7 inventory can be 
found. 

Cumulative mass depth The distribution of 7Be with depth in the soil 


profile is expressed in units of areal mass 
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(kg m7”) rather than metric units (mm). This 
provides a more accurate measure of depth in 
the context of the soil redistribution models. 

Profile Distribution Model (PDM) Simple model to convert 7Be inventory loss 
into soil erosion amounts by linking inventory 
change to the depth profile to derive the mass 
depth of soil loss. This model is applied to sin- 
gle events over a short period of time. 


Appendix 1.1: The Profile Distribution Model 


The exponential depth distribution of Be in a soil profile can be described as: 
C(x) = ce*/"0 (1.1) 


where C, (Bq kg~') is the ’Be activity at mass depth x (kg m~”), c is a constant 
value and hg is the relaxation mass depth, the depth at above which 63.2% of the ’Be 
inventory can be found. 

Erosion rates (kg m~7) can be estimated by comparing the ’Be inventories at the 
sample site, A (Bq m~7), to the reference inventory, Aref (BQ m~*). Where a mass 
of soil has been lost (h) (kg m~?) changes in the sample site inventories can be 
represented as: 


CO 
A= f C(x)dx = Arefe” ™ (1.2) 


—h 
Erosion rate, h (kg m~, negative), can, therefore, be calculated as: 
h = ho ln(A/Aref) (1.3) 


Deposition of material is reflected in an excess of ’Be inventory at the sample site 
with respect to the reference site. The depth of deposition, h’ (kg m~’, positive), can 
be calculated as: 


h' = (A — Aver) /Ca (1.4) 


Where C4 (Bq kg!) is the mean activity concentration of ’Be in the deposited 
sediment. 

The ’Be activity concentration in the eroding sediment at each upslope point, C, 
(Bq kg~'), can be calculated from the loss of inventory divided by the mass of soil 
loss: 
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Ce = (Aref — A) /h (1.5) 


The mean activity concentration of soil mobilised from the study area, S (m7), 
can then be calculated as: 


ca= f nc.asy f nas (1.6) 
S 


S 
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Chapter 2 A) 
How to Design a Be-7 Based Soil get 
Distribution Study at the Field Scale: 

A Step-by-Step Approach 


W. H. Blake, A. Taylor, A. Toloza and L. Mabit 


2.1 Key Sample Sets and Associated Data 


A 'Be-based soil redistribution budget is based on several key datasets (Table 2.1) 
which have strict rules on collection locations and timings, depending on the time 
period of application. Overall, the methodological approach follows the principles 
of other FRN techniques (e.g. the '°’Cs approach) but with necessary differences 
linked to the short half-life of 7Be and its delivery dynamics. The difference in 
delivery dynamics also provides the added advantage of opportunity for assumptions 
underpinning the approach to be tested in field and by laboratory experimentation 
(Taylor et al. 2014). As summarized in Table 2.1, some datasets are mandatory 
to convert measurement of Be inventory into soil redistribution amounts. Other 
datasets are advised under some circumstances to assist with data interpretation and 
improve the quality of soil distribution estimates. Fundamental considerations for 
the collection of all these datasets are outlined in the following section. 
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2.2 Reference Site Selection and Sampling 


As outlined in Chap. 1, ’Be inventory data are used to construct an FRN budget 
for the hillslope, to evaluate relative differences to the reference inventory, and this 
budget is subsequently converted into a soil redistribution budget (described in detail 
in Chap. 4). Accurate and representative determination of the reference inventory, 
and associated uncertainty is therefore a fundamental requirement since estimates of 
soil redistribution pattern and amount relies on this key value. 

Reference sites serve two purposes in ’Be studies: (1) to determine a mean refer- 
ence inventory (i.e. the areal activity of "Be (Bq m~”) in the soil surface unaffected 
by erosion), and (2) to determine the depth distribution of Be in the undisturbed soil 
prior to erosion, to derive ho. The first requirement is based on collection of bulk 
soil cores, and the second on the collection of at least one sectioned core. Both need 
to be determined at a flat, stable location near the study plot, as for '°’Cs studies, 
but it is essential that the land-use history of the location for the depth profile is 
exactly the same as the study plot. This is important as the ho value determined must 
be representative of the eroding soil surface. In practice, the best location for both 
these measures is a flat area that has been cultivated at the same time and in the same 
manner as the study plot (Fig. 2.1). 

Sample designs for bulk cores need to account for potential FRN spatial variability 
within the reference site (Sutherland 1996; Mabit et al. 2012; Kirchner 2013) and 
also deliver sufficient mass of soil for sample analysis (Chap. 3). However, sample 
numbers are often constrained in Be tracing investigations owing to the short half-life 
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Fig. 2.1 Experimental design to establish a hillslope 7Be budget over largely bare soil 
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and availability of sufficient gamma detection facilities. Considering this limitation, 
Taylor et al. (2013) recommend that all studies should state a reference inventory 
with suitable upper and lower limits (i.e. +20) (Owens and Walling 1996) and this 
should be incorporated into subsequent soil redistribution estimates (Chap. 4). Spatial 
variability in inventory within a plot is most likely due to local redistribution of soil 
by rain splash and micro-topography and accumulation of rain-splashed particles in 
hollows. To capture such variability at any given reference site, it is recommended 
that each sample integrates several cores and then this process is replicated to provide 
a minimum of 10-15 spatially-integrated reference samples for analysis (Fig. 2.1). 
Bulk soil cores need to be taken to a consistent depth that is below the known 
depth penetration of ’Be in the study soil. Depth penetration is typically 20-30 mm 
(Blake et al. 1999; Doering et al. 2006; Sepulveda et al. 2008; Wallbrink and Murray 
1996). In this regard, we strongly advise undertaking a preliminary investigation for 
trial depth-profile dataset at the study area to ensure the complete profile is captured 
whilst avoiding dilution of the sample by overestimating the profile depth. This will 
allow the researcher (1) to establish the maximum depth penetration of Be in the soil 
of the investigated area and (2) to use this depth to determine the maximum depth 
collection of the remaining bulk cores to be collected (e.g. 0-30 mm in the reference 
site). Any vegetation on the ground surface must be included in the sample as this will 
carry part of the recent ’Be inventory (Iurian et al. 2015). When collecting spatially- 
integrated cores, it is important that the total sum of the core areas that comprise one 
sample is recorded and that the depth penetration of all cores is consistent at 1 mm 
precision. 

When characterising the ’Be depth profile, it is essential that section cores are 
sampled from a soil surface that has experienced the same cultivation practice as the 
study slope but in a location that has remained undisturbed by erosion or deposition 
processes. Variability due to rain splash, as noted above, can present challenges in 
selection of the appropriate position of the core. It is highly recommended that more 
than one section core is collected but this is often limited by gamma detector resource 
availability. An alternative is to collect 3 replicate section cores and combine the 
respective layers from each core into integrated samples to capture, but not to quantify, 
spatial variability within the flat, uneroded reference area. The recommended tool 
for standardized section core sample collection is the Fine Increment Soil Collector 
(FISC) (Mabit et al. 2014; Fig. 2.2) which is proven to collect high precision depth 
profile data suitable for supporting "Be inventory conversion (Ryken et al. 2016). 

It is also common practice to establish a rainfall collection and monitoring station 
at the reference site to permit assessment of the dynamic of ’Be delivery in relation to 
inventory development and its radioactive decay in the study area under investigation 
(Fig. 2.3). Sampling of reference cores through a time period also serves to benchmark 
and/or validate rainfall-based inventory assessment (Wallbrink and Murray 1994; 
Walling et al. 2009). The importance of this when applying the event scale PDM is 
to validate the stability of the reference site. With high resolution rainfall data and 
rainwater samples, we can use the modelled inventory as a benchmark for confidence 
in our choice of reference site. For this purpose, rainfall data need to be collected 
using a tipping bucket rain gauge (Fig. 2.4) that provides 15 or 30 min interval 
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Fig. 2.2 Fine Increment Soil Collector (Mabit et al. 2014) 


rainfall intensity data. Alongside rainfall monitoring, bulk rainfall samples need to 
be collected at a minimum of monthly intervals but preferably at event-scale intervals 
(depending on rainfall regime) and analysed for "Be concentration (Bq L~'). Samples 
should be collected in pre-acidified polypropylene containers with an attached funnel 
for rainfall capture (Fig. 2.3) and Be extracted following protocols described by 
(Taylor et al. 2016); (see Appendix 2.1 for standard operating procedure). 

A time series of the relative inventory at a study site (e.g. Fig. 2.3) can then be 
calculated as follows, where for each daily time step: 


A(t) = A(t — 1) - exp(—A) + F(t) (2.1) 


where A(t — /) is the Be inventory of the previous day (Bq m7”), » is the radioactive 
decay constant (daily time step) and F(t) is the fallout contribution of the current day. 
It is common practice to determine the start point either through collection of a suite 
of reference inventory cores (Walling et al. 2009) or to assume zero due to tillage of 
the soil surface and dilution of the "Be signal within the soil profile. Table 2.2 shows 
a spreadsheet coding example to create a cumulative inventory dataset. 
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Ī Daily rainfall 


— Be inventory 


7Be inventory (Bq m?) 
Daily rainfall (mm) 


15/11/97 01/12/97 15/12/97 01/01/98 


Fig. 2.3 Example of reference inventory development and decay record based on rainfall sampling 
and analysis for "Be alongside rainfall monitoring (Blake et al. 1999) 


Ses 


Fig. 2.4 Examples of tipping bucket rain gauge (left) and rainfall sampler (right) 


2.3 Sample Design Options for Soil Redistribution 


The ’Be approach is most suited to erosion studies performed on bare soil surfaces 
(no vegetation) or with little vegetation cover. The method is limited to quantification 
of erosion by rain splash, sheet wash and shallow rill development since once rill 
incision goes beyond the depth of the "Be depth profile, eroded soil is exported in the 
absence of the tracer signal. The selection of study site is very much dependent on the 
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Table 2.2 Coding for Microsoft Excel spreadsheet to calculate cumulative inventory from rainfall 
data and rainwater ’Be activity concentration data 


Cell number (assuming first Code Explanation 
row is column headings NB 
there is only one entry on row 
2—see below) 


A3 n/a Date and time step of rainfall 
data 

B3 n/a Rainfall input for time step 
(mm) 

C3 =B3 x 1000/1000 Rainfall expressed as volume 
per unit area (L m~*) 

D3 n/a Rainfall 7Be activity 


concentration (Bq L7!) 


E3 D3 * C3 7Be areal activity deposition 
for time step (Bq m7?) 


F2 n/a Inventory at start of 
monitoring (generally zero if 
plot is cultivated at beginning) 


F3 =E3 + (F2 x EXP(-$I$1)) Plot inventory (Bq m7”) at 
timestep where cell I1 
contains the decay constant 
for the appropriate time step. 
The ’Be deposition received 
(cell E3) is added to the decay 
corrected inventory from the 
previous day 


research question but, in any case, a basic pilot study will serve useful in determining 
the presence of a sufficient ’Be inventory in the region to make the approach viable 
and, as described above, to evaluate the typical depth penetration of ’Be to inform 
bulk core sample depth. 

The "Be approach may be applied at different spatial scales depending upon the 
questions being asked by land managers and the constraints imposed by the key 
assumptions discussed in Chap. 1. In this context, there are three main ways to 
design a soil core sampling strategy within the study hillslope: (1) sampling along 
transects from upper to lower slope (Schuller et al. 2006); (2) high resolution grid 
sampling (Walling et al. 1999; Blake et al. 1999); (3) spatially-integrated sampling 
within defined geomorphic landscape units (Wallbrink et al. 2002; Blake et al. 2009). 

At the field scale, the single transect (Fig. 2.1) or a multi-transect sampling is the 
most straightforward and cost-effective approach in terms of field sample collection 
and laboratory processing work effort. Transects can, however, be limited in terms of 
spatial representativeness depending on local topography and research questions to 
be addressed. A high-resolution regular grid approach is effective for evaluation of 
spatial variability and will provide more representative information on soil export and 
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sediment delivery ratio from a larger study area than transects. However, this approach 
is highly limited by the analytical demand of a large number of samples, which 
is hampered further by the short half-life of 7Be. The geomorphic landscape unit 
approach offers a compromise between the limited spatial representativeness of the 
transect approach and the analytical demands of the grid approach. However, it should 
be pointed out that multiple reference sites may be required in some geomorphic 
landscapes. 


2.4 Sampling for Particle Size Selectivity Correction 


Particle size selectivity during soil erosion processes is well-known (Bernard et al. 
1992) and users of soil loss assessment or measurement approaches need to make 
a decision on the relevance of this specific process to their own study site. Taylor 
et al. (2014) describe how particle size selectively of soil erosion processes can be 
accounted for in the 7Be conversion model. Application of this method requires 
specific samples to be collected during or after the erosion event being studied. 

Particle size correction requires representative samples of (1) uneroded soil, (2) 
mobilised soil and (3) deposited soil to be collected (Fig. 2.5). The first can simply 
be represented by the bulk cores collected for depth profile determination (i.e. in 
a non-eroding but cultivated site). The third requirement above can be represented 
by the samples collected to determine inventory in areas of sediment accumulation. 
The second, the characterisation of mobilised soil, requires more careful planning 
since after an erosion event, such material has either left the study site or been, in 
part, deposited at the foot slope. Taylor et al. (2014) propose that such material can 
be collected with Gerlach troughs installed in the study site or simple equivalent 
sediment trap systems. Alternatively, to avoid the need to install equipment prior to 
the event, a rainfall simulator could be used to mobilise material from a representative 
area, and the material captured for analysis. This needs to be done under rainfall 
conditions similar to the erosive event. 

All samples collected need to be analysed for particle size distribution (Chap. 4). 


2.5 Summary: Designing a Basic Small Scale ’Be Pilot 
Study Sampling Programme at the Plot Scale 


The procedures described in this chapter can be practised and developed within the 
context of a simple pilot study framework designed to test the viability of using “Be 
as a tracer in any given landscape. A step-by-step approach is detailed below: 


Step 1: Having located your study site and secured permission from the landowners 
for its investigation, collect 3 topsoil samples (20 mm depth) after significant rainfall 
(e.g. >sufficient rainfall to develop a measurable inventory for the study area over an 
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extended time period) and analyse them for ’Be content to confirm if significant and 
measurable inventory of this radioisotope is present in the study area; 

Step 2: Collect a trial section core to evaluate the depth penetration of "Be in the 
study soil; 

Step 3: Select the reference site in the study area and, if possible, the study plot 
location (prior to soil erosion taking place). Set up a rain gauge and, if required, rain 
sampling equipment. Note the timing of the last cultivation and, if this was some time 
prior to current time within which rainfall had occurred, collect a set of reference 
cores to establish the baseline (t = 0) inventory. Note this will be zero by default if 
plot study begins immediately after cultivation and mixing of existing ’Be inventory 
into the soil profile; 

Step 4: After erosion has taken place within a target study plot (a) collect 10-15 
spatially-integrated bulk reference cores from the non-eroding reference site to the 
depth of maximum ‘Be penetration (e.g. 20-30 mm) following the information pro- 
vided by the test core (in step 2), (b) collect up to 3 section cores from a flat, 
non-eroding area that has experienced the same cultivation practice as the study 
plot using the FISC and either analyse cores separately or combine layers to cre- 
ate a spatially integrated depth profile depending on analytical resource, (c) collect 
spatially-integrated cores along three transects within the study plot (Fig. 2.1) with 
5-15 samples per transect depending on your analytical resource; 

Step 5: Undertake necessary sampling for particle size correction if desired; 

Step 6: Bring all recovered soil samples to laboratory (i) for preparation prior to 
gamma spectrometry (Chap. 3) measurements and, when necessary, (ii) for perform- 
ing particle size distribution analysis by laser granulometry. 


* Fine 
aaa + Medium 


Coarse 


e 

S Uneroded soil 
m Mobilised soil 
d Deposited soil 


Fig. 2.5 Schematic diagram of sampling protocol for including particle size selectivity in ’Be soil 
erosion study (uneroded soil = s; mobilised soil = m; deposited soil = d) 
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Glossary 

Reference Inventory Radioactivity per unit area (Bq m7”), also termed 
areal activity, at a stable (non-eroding) field site in 
close proximity to the eroding study site. 

Bulk soil core A soil core taken to a specific depth wherein the 
material recovered within the corer represents one 
sample. 

Section core A soil core that is initially kept intact within the soil 


corer and then extruded in small (ca 2 mm) incre- 
ments to permit sections to be subsampled layer by 
layer. 

Spatially—integrated sampling A process through which one sample collected for 
analysis comprises several smaller samples taken 
over a wider area to capture local variability e.g. due 
to micro-topography or variable vegetation cover. 
Can be used in routine core sampling but also 
extended to underpin the landscape-unit sampling 
approach to "Be budgeting (see text for details). 

Particle size selectivity The process by which erosion processes prefer- 
entially remove fine-grained soil particles due to 
greater critical shear stress required for mobilisation 
of coarser particles. 


Appendix 2.1: Protocol for Extraction of "Be from Rainwater 


According to Taylor et al. (2016), prior to deployment of rainfall collectors, 10 mL 
HCI (2.5 M) is added to each bottle to prevent adsorption of ’Be to vessel walls during 
the sampling period. Due care must be taken when handling acid in accordance with 
your institutions Health and Safety code. Collectors may be exposed for periods 
of between 3 and 35 days depending on the frequency and magnitude of rainfall. 
Generally, samples comprise fallout from a number of events across the sample 
period and are therefore, referred to as integrated samples. At the point of sampling, 
funnels should be rinsed with a known volume of HCI (1 M) and the bottles replaced 
with acid-cleaned bottles. 

Each sample should be checked to ensure pH is < 2 and then filtered to remove any 
coarse debris (e.g. using Whatman grade number 41 filter papers). 7Be can then be 
pre-concentrated from solution by co-precipitation with MnO; following the method 
detailed by Short et al. (2007). 

1 mL of 0.2 M KMnO; is added per litre of rainwater sample and the pH adjusted 
to 8-10 using concentrated NH4,OH. 
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Once at the desired pH, 1 mL of 0.3 M MnCl, is added to the sample whilst 
stirring. MnO% precipitate is then allowed to settle for 24 h prior to removal by 
vacuum filtration using 0.45 um cellulose nitrate filter paper. 

Filter paper is then air dried, fixed with cellophane and sealed in a suitable con- 
tainer prior to analysis by gamma spectrometry. For each rainwater sample, duplicate 
1 L subsamples should ideally be treated and the precipitate combined for filtration. 
Where rainfall samples are of low volume, a single 1 L sample may be treated. Sam- 
ples should be considered to provide a representative sub-sample of rainfall for the 
period. 

Repeatability can be tested by analysing triplicate subsamples separately and Rel- 
ative Standard Deviation (RSD) between triplicates determined. This should typically 
be <10% (Taylor et al. 2016). 

7Be recovery from solution using the coprecipitation method was tested by Taylor 
et al. (2016) by reprecipitating the filtrate from 3 samples. In each case, the ’Be 
activity in the filtrate was below Minimum Detectable Activity (MDA). MDA values 
for these samples were <10% of the total activity, suggesting that ’Be recovery using 
the above method is >90%, in agreement with Short et al. (2007). 
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Chapter 3 A) 
Measurement of Be-7 in Environmental E 
Materials 


A. R. Iurian and G. E. Millward 


3.1 Overview of Digital Gamma Spectrometry Systems 


Gamma-ray spectrometry is the only ‘routine’ method for measuring the natural cos- 
mogenic radionuclide "Be. Activity concentrations (Bq kg~') of ’Be are determined 
by analysing the 477.6 keV gamma energy emitted by the first excited state of 7Li 
(branching ratio 10.44%) (DDEP 2017) as it achieves the ground state. A typical 
gamma-spectrometric system consists of a semiconductor crystal detector, liquid 
nitrogen or mechanical cooling system, preamplifier, detector bias supply, linear 
amplifier, analogue-to-digital converter (ADC), multi-channel analyser (MCA) of 
the spectrum, and output data devices (ANSI 1999). An example of a typical gamma 
detector is given in Fig. 3.1. 

High purity germanium (HPGe) detectors currently represent the most widely used 
gamma-detector systems. The semiconductor crystal is manufactured from ultrapure 
germanium (impurity level 109 atoms cm~*) in various shapes and in a range of sizes 
for a wide range of applications (Debertin and Helmer 1988). The specification of 
a HPGe detector is defined by (i) the energy resolution, (ii) the detection efficiency 
and (iii) the peak-to-Compton ratio (see Glossary). When purchasing a HPGe semi- 
conductor detector, these key parameters need to be considered together with the 
detector material and configuration, detector volume, window material and its thick- 
ness. Excepting the ultra-low energy configurations, all detector types commercially 
available are suitable for the measurement of ‘Be activity concentrations. 
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(b) 


Fig. 3.1 a Example of a HPGe gamma detector showing the detector head, the cartridge containing 
electronic management system and the dip-stick which is immersed in liquid nitrogen to cool down 
the detector to 77 K; b the detector encased in a lead shield with the dip stick immersed in a liquid 
nitrogen tank. NB: The liquid nitrogen tank should be filled on a regular basis (every 7-10 days) 
to maintain a fully operational detector 


3.2 Preparation of Samples 


3.2.1 Sample Management 


A fundamental consideration for the analyst when managing a set of soil, sediment 
or rainwater samples is the relatively short half-life of Be and its low activity in the 
majority of environmental materials. It is important, therefore, that the laboratory 
is well organised and ready to process the samples for 7Be analysis as soon as 
possible after collection so that a relatively high emission of gamma rays can be taken 
advantage of. Additionally, maintenance of a high standard of cleanliness in order 
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to prevent contamination of the detector end cap or potential inter-contamination 
between samples should be prioritized. 


3.2.2 Soil Samples 


Samples of soil will normally consist of a range of particle sizes. Gamma spectromet- 
ric analysis of fallout radionuclides (FRN) in bulk soils should typically be performed 
on dried, homogenised material that has been sieved to less than 2 mm (Pennock and 
Appleby 2002) noting that some applications may require sieving to a smaller par- 
ticle size. Typically, soil samples are oven-dried at 105 °C until a constant weight is 
reached (commonly for 24 h). A practical alternative is sample freeze-drying, which 
disaggregates the solids and makes subsequent sample processing more effective. 
Furthermore, oven-dried samples are disaggregated using a mortar and pestle and 
size fractionation of the dried solid can then be carried out on the ground sample 
using a metal or plastic sieve with a defined mesh size. Appendix 3.1 summarises 
the sample preparation and the associated packing protocol. 


3.2.3 Samples of Rainwater 


Samples of rainwater from the study site are required to assess the atmospheric 
deposition flux of "Be during a period, encompassing the sampling campaign (see 
Chap. 2). The protocol for extraction of ’Be from rainwater is presented in Appendix 
2.1. 7Be can be pre-concentrated from solution by co-precipitation with MnO 
following the method detailed by Short et al. (2007). The precipitate is recovered on 
glass fibre filters which are placed in a Petri dish of suitable diameter and allowed 
to dry before counting on an HPGe gamma detector. 


3.2.4 Sample Geometry 


Subsequent to physical processing samples, for analysis, can have masses in the range 
from a few grams to several kilograms, according to the aims of the project and its 
sample collection requirements. There are a range of measurement geometries, from 
Marinelli beakers (0.25—1 L) to small vials (~5 cm?) capable of holding a wide range 
of solid weights (from about 3 g of material to 4 kg) (Fig. 3.2). The choice of the 
container geometry for "Be analysis depends on the matrix and on the available 
amount of material. It is essential that the weights of both empty and filled sample 
holders are registered on a calibrated laboratory balance to at least 2 decimal places 
(or at least 4 decimal places for volume calibration standards). 
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Fig. 3.2 Several examples of containers used to perform gamma spectrometric analysis of solid 
phases. From left to right: 1 L Marinelli beaker; 300 mL Marinelli beaker; 90 mm diameter plastic 
Petri dish; 50 mm diameter Petri dish; 4 mL plastic vial (commonly used for conducting analysis 
using well detectors). To comply with laboratory safety, the radioactive standards are bound in red 
tape 


Prior to purchasing any sample/standard container (especially for Marinelli beaker 
geometries), the analyst should determine the diameter of the detector end cap so that 
the container with an appropriate internal diameter will fit as close as possible on the 
detector end cap. For cylindrical and vial containers, it is recommended to place the 
samples on a holder (e.g. of light polycarbonate matrix) to assure the reproducibility 
of the sample-detector arrangement for all samples, including the measurement of 
the standard. The use of such holders will prevent contamination of the detector end 
cap, which can be further protected with a thin sheet of cling film or, in the case of 
well detectors, the sample vial covered in film. Care should be also taken to ensure 
that the sample/standard container external surface is free of particulate matter. 

It is recommended that the geometries of the standard and the soil samples are 
identical; any differences between them, such as in the sample filling height (e.g. 
because of insufficient sample amount) or in soil chemical composition, will imply 
the need for geometry corrections (see Glossary). It is important to note that if the 
containers are being re-used, the analyst needs to carry out acid cleaning and drying 
procedures before refilling them. The packing date, nominal weight and an identifi- 
cation code should be indicated on each sample container after packing. As regards 
filter-mounted samples, these should be measured after placing them in cylindrical 
Petri dishes of appropriate diameter. Similar contamination precautions are to be 
taken for the filter geometry as for the other container types discussed above. 
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3.3 Calibration Approaches 


An accurate and precise calibration of a gamma detector is essential for high qual- 
ity results in the determination of ’Be activity concentrations. The calibration of a 
gamma spectrometric system involves two major steps: 


(1) Energy calibration (includes energy calculation as a function of the number of 
the channel and Full-Width-Half-Maximum (FWHM) as a function of energy); 

(2) Efficiency calibration (the efficiency as a function of the energy correlated with 
the source geometry). 


Energy calibration is commonly performed using point sources (e.g. !°*Eu, !°’Cs) 
or multinuclide standards covering the energy range of interest (Fig. 3.3). The con- 
struction of the full energy peak efficiency curve is performed following (1) a 
semi-empirical approach (based on measurement of standard point sources or multi- 
nuclide standard sources of soil or liquid matrix), or (2) Monte Carlo methods (using 
a detailed description of the detector and the source). In both cases, the efficiency 
calibration procedure involves (1) true photopeak efficiency determination for the 
gamma energies of the radionuclides included in the certified standard source, and 
(2) the construction of the efficiency curve (within a range of energies given by the 
radionuclides present in the source) by polynomial fitting. Discrepancies between 
simulated (Monte Carlo) and experimental efficiency values should be determined 
and included as an uncertainty component in the uncertainty budget of the ’Be activ- 
ity value. Note that different detectors and different system-sample configurations 
will result in different values of the true photopeak efficiencies, which should finally 
produce the same result of the radionuclide activity value for the actual sample. 
Coincidence summing corrections should be performed for relevant nuclides present 
in the calibration source and decaying through a cascade of successive photon emis- 
sions (e.g. ®Co, 88Y, !3°Ce). These can be determined through general or dedicated 
codes (Vidmar et al. 2010). Coincidence summing correction factors depend on the 
nuclide decay scheme, sample geometry and composition, and on detector parame- 
ters. The factor is equal to unity in cases where the radionuclide has no cascade of 
gamma-rays. 


3.3.1 Semi-empirical Calibration Method 


For routine environmental sample analyses, efficiency calibration of gamma detec- 
tors can be performed using (a) multiple standard point sources or standard 
volume sources of mixed radionuclides (for soils measurements), or (b) filters 
spiked with standard solution (for 7Be precipitate analysis). Considering the 
short half-life of 7Be, standard sources containing "Be are not commercially 
available. When purchasing a standard for efficiency calibration, it is essential 
that these should carry a certificate of calibration from a metrological labora- 
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Fig. 3.3 Print screen capture of the gamma detector energy calibration performed using a multi- 
radionuclide standard and the GammaVision commercial software. Note the interpolated 477.6 keV 
gamma energy of "Be between 391.7 keV (!!3Sn) and 661.7 (!37Cs). The upper box is the energy 
table; the middle box is a plot of the gamma-ray energy against the channel number; the lower box 
shows the 7Be peak 


tory accredited for the preparation of standards. An appropriate calibration stan- 
dard should have low overall uncertainties of each radionuclide in the source 
(up to 2~-3%), as this uncertainty is further propagated into the uncertainty of ’Be 
activity. The analyst should be aware that specific international regulations apply for 
the import of radioactive sources. After purchasing, the calibration standard solution 
(commonly delivered in sealed glass vial) must be gravimetrically diluted for sub- 
sequent use. Furthermore, the analyst shall choose how to use the diluted standard 
solution for the detector efficiency calibration. Bear in mind that special safety pro- 
cedures and authorisations apply to the laboratories handling radioactive sources as 
regards to the risk of radioactive contamination. Appendix 3.2 details arecommended 
step-by-step approach for the preparation of a multi-nuclide calibration standard of 
soil matrix. 

An example of a multi-nuclide standard solution suitable for Be measurements 
will comprise artificial radionuclides such as **! Am, !°Cd, *’Co, !3°Ce, ?Hg, !Sn, 
85Sr, 137Cs, Co and 8Y, with energy lines from 59 to 1836 keV. Such multi-nuclide 
standard can be used to calibrate gamma detectors for efficiency in the vicinity of the 
gamma emission from Be (477.6 keV) which lies between that of !!?Sn at 391.7 keV, 
85Sr at 514.0 keV (difficult to determine accurately due to the overlapping with the 
511 keV annihilation peak) and !3’Cs at 661.7 keV (DDEP 2017), as presented in 
Fig. 3.4. A disadvantage in the short half-lives of some nuclides (e.g. 7°°Hg, !'Sn) 
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Fig. 3.4 The efficiency calibration curve of a HPGe detector using a multi-nuclide standard solution 
containing 10 artificial radionuclides (plotted in Microsoft Excel spreadsheet). Note the 7Be gamma 
emission is interpolated between the gamma energies of '!3Sn and ®*Sr. After the decay of short- 
lived radionuclides i.e. Sr, !!3Sn and 2°Hg, the closest lower energy for interpolation will be 
'39Ce, An interpolation of "Be energy between !3°Ce and !37Cs will led to an increased efficiency 
value for 7Be 


present in the standard must be recognized, thus making the source useful over the 
whole spectra region only for a short period of time (i.e. a few months). 


(i) Preparation of a secondary soil calibration standard 


A soil sample should be collected, preferably from an area where fallout radionu- 
clides are minimal (e.g. subsoil) and prepared as in Sect. 3.2.2. The detection and 
quantification of the background gamma radiation in the soil sample is required. 
Then the prepared soil fraction below 2 mm particle size is spiked with an aliquot 
of the liquid standard. It is highly important to ensure that the soil-liquid standard 
mixing is complete, ensuring the spiked material is homogeneous (see step-by-step 
instruction in Appendix 3.2). The spiked soil is dried and an appropriate sample 
holder (selected from those in Fig. 3.2), with a similar geometry to the samples, 
is then selected and filled with the dried spiked soil. The secondary soil standard 
can then be used to calibrate the detector for efficiency (Fig. 3.5). If there are dif- 
ferences in geometry (different containers and/or filling height) and/or matrix (e.g. 
different chemical composition) between the soil samples and the calibration source, 
an efficiency transfer factor should be determined to correct this difference, as ratio 
between the efficiency of the sample and the efficiency of the standard. 


(ii) Preparation of a secondary liquid calibration standard 
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Fig. 3.5 Print screen capture of the gamma detector efficiency calibration performed using a multi- 
nuclide standard and the GammaVision commercial software. The upper box is the energy table; 
the middle box is a plot of the efficiencies against the gamma-ray energies; the lower box shows 
the Be peak 


In this case, the liquid standard commonly purchased in a glass vial is diluted follow- 
ing the steps in Appendix 3.2, and further poured in different geometries and used 
as a secondary liquid standard. In this case, the efficiency transfer factor needs to be 
considered for ‘Be activity calculation in soil samples, to account for the difference 
in matrix between the calibration source and the sample (regardless the sample filling 
height). 


Gii) Preparation of a secondary filter standard 


A secondary filter standard can be obtained by spiking a filter paper in similar size 
with the sample (e.g. Whatman filter) with a known content standard solution. This is 
achieved using a calibrated automatic pipette that is used to add the standard solution 
dropwise onto the filter. Here the analyst must take care to ensure a homogenous 
distribution of the standard solution on the filter so that the calibration is coherent 
with the distribution of precipitate on the sample filter. The spiked filter should then 
be air-dried prior to sealing in a plastic Petri dish of appropriate diameter. 


3.3.2 Monte Carlo Approach 


The advanced user may prefer to use Monte-Carlo coding to construct the efficiency 
curve for a certain sample holder without the need of a certified standard solution. 
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Monte Carlo approaches for automatic efficiency calibration (e.g. LabSOCS—Lab- 
oratory SOurceless Calibration Software) (Bronson and Wang 1996) are commonly 
integrated into some commercial software of gamma detectors. These codes include 
several options for detector type and configuration, sample shape and volume (e.g. 
point, disc, and cylinder), sample/source matrix, source-to-detector distance and for 
the type of fit used to describe the efficiency-energy dependency (linear, quadratic, 
polynomial) (Jovanovic et al. 1997; Bronson and Wang 1996). For each particular 
sample, the analyst needs to create a specific sample geometry based on the holder 
dimensions (diameter and wall thickness), sample filling height, chemical composi- 
tion of the sample matrix and density. It is also important for the analyst to know the 
precise details of the technical characteristics of the gamma detector or to obtain a 
‘characterisation file’ of the detector from its manufacturer in order to run the codes. 
However, no programming knowledge is necessary to use the available commercial 
codes. The efficiencies generated for a specific sample can be saved and stored by 
the analysis software and used to automatically derive the radionuclide activity with- 
out performing an experimental calibration. However, the experimentally obtained 
efficiencies can be as well included in the analysis software, in case these data are 
available. 

The automatic calibration by the Monte Carlo approach eliminates the cost of 
purchasing, tracking, and disposing of radioactive standards. However, the analytical 
results still need to be validated using certified reference materials in similar matrices. 
The main constraints in the application of Monte Carlo approaches are related to 
uncertainties of the detector configuration (e.g. shape and size of the effective crystal 
volume, photons and electron interaction probability and angular distributions) and 
measurement traceability to primary standards. 


3.4 Data Handling and ’Be Analysis 


Analyses of ’Be can be challenging because of low activity concentrations and, 
therefore, high counting uncertainties. Sample measurements are completed after 
performing the energy and efficiency calibrations. Data evaluation is commonly 
realised using specialised gamma software, provided by the detector manufacturer 
or made ‘in-house’. Prior to spectrum analysis, 7Be nuclear data (gamma energy, 
half-life and emission probability) should be included in the software library. The 
stability of the spectrometer should be checked during the measurement and during 
routine quality assurance procedures. Despite the apparent simplicity of gamma- 
ray measurements, bear in mind that there are a number of correction factors to 
the spectrum counting data that must be considered for the determination of ’Be 
activity concentration e.g. decay correction to the sampling date, efficiency transfer 
factor, background correction, instrument dead time correction IAEA 2004), while 
a correction for soil moisture content is applied to the dried soil sample mass. 

A flow chart of 7Be measurement by gamma-ray spectrometry is presented in 
Fig. 3.6. Sample counting times need to be long enough (typically between 86,000 
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Fig. 3.6 Schematic protocol for 7Be analysis by gamma-ray spectrometry (MDA = Minimum 
Detectable Activity) 


and 259,000 s) to meet laboratory standard statistical uncertainty of the measurement 
result (commonly between 10 and 30% for ’Be in environmental samples). 

Based on the information provided in the library, the instrument-specific gamma 
evaluation software will identify the 7Be peak and analyse the number of counts 
within its peak area (Fig. 3.7). Furthermore, following default automated sequences; 
the software is able to report the activity concentration (Bq kg~!), combined standard 
uncertainty and minimum detectable activity of the radionuclide found by the peak- 
search routine (Debertin and Helmer 1988). However, to achieve reliable results 
using the library approach, any shift effect in the spectrum should be corrected and 
energy calibration should be precise enough to achieve no more than one or two 
channels deviation between ’Be reference energy and the one calculated from the 
calibration curve (IAEA 2002). 
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Fig. 3.7 Print screen capture of a gamma spectrum registered with the GammaVision commercial 
software and showing the ’Be peak at 477.56 keV and its net area and gross/net count rate (peak 
data box), measurement starting time and live time (upper right part of the screen) 


Measurement results are commonly reported with associated combined uncer- 
tainty at the 95% confidence level (JCGM 2008). For detailed information on the 
determination of each uncertainty component, readers can refer to Lépy et al. (2015). 
The detection limit (LD) provides the detection capabilities of a measurement system. 


3.4.1 Quality Assurance 


In routine laboratory practice, it is necessary to validate the measurement results 
for certain working conditions. Both internal and external validation methods are 
an important part of ‘good laboratory practice’ and also a requirement of ISO/IEC 
17025 for ‘in-house’ methods (ISO/IEC 17043 2010; ISO/IEC 17025 2017). The 
internal validation procedure can be employed using distinctive approaches: 


— Various types of blanks; 

— Replicate analysis to check changes in precision; 

— Use of certified reference materials with similar matrix and geometry (e.g. provided 
by well-recognized international providers such as the IAEA); 

— Standards and point sources to check the stability of the response; 

— Use of long-term control charts of standards and point source measurements to 
check the stability of the detector response; 

— Proficiency tests participation. 
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Note that because of its short half-life, there are no reference materials available 
which contain ’Be. However, it is a common practice to check the accuracy of the 
laboratory procedure for radionuclides emitting gamma rays in the vicinity of Be 
energy (e.g. !37Cs) in samples of similar matrix and geometry. 

In addition, external validation can be achieved through: (i) inter-laboratory 
comparison exercises and/or (ii) participation in international proficiency tests for 
radionuclide determination in soil and water samples such as the ones organized 
regularly by IAEA and its Member States. 


Glossary 


The energy resolution (FWHM) The full width of the peak at one- 
half of the maximum height above 
any underlying continuous spec- 
tral background.! 

Absolute efficiency The ratio between the number of 
pulses recorded and the number 
of radiation quanta emitted by 


source. 

Relative efficiency Detector efficiency relative to that 
of a Nal(T1) detector.” 

Peak-to-Compton ratio The ratio of the number of counts 


in the biggest channel of the 
1332.50 keV Co energy peak 
to the average channel count in 
the Compton continuum between 
1040 and 1096 keV in the same 
spectrum. 

Coincidence summing correction factor The ratio between the appar- 
ent efficiency for the line with 
energy E of the nuclide hav- 
ing cascading radiations, to the 
full energy peak efficiency at 
the same energy obtained from 
the energy curve measured with 
single-photon emitting radionu- 
clides. The apparent activity (cor- 
rected for background) must be 
multiplied by the coincidence 
summing correction factor to 


'Debertin and Helmer (1988). 
Knoll (2000). 
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Geometry and self-attenuation correction factor 


Efficiency transfer 


Minimum Detectable Activity (MDA) 


Detection Limit (LD) 
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obtain the true activity value for 
the radionuclides affected by cas- 
cade effect.* 

The ratio of the full energy 
peak efficiency for the sample 
and the full energy peak effi- 
ciency for the calibration stan- 
dard. The correction is commonly 
performed through an efficiency 
transfer approach, by multiplying 
the correction factor with the ’Be 
full energy peak efficiency. 

The “transfer” of the efficiency 
of a standard to another sample, 
which can be in the same geome- 
try or not, and acting as corrective 
factor for both geometry and self- 
absorption effect. 

The lowest activity value that can 
be detected with a measuring sys- 
tem at 95% confidence.* 

The lowest expectation value of 
the counting rate of a net peak area 
that can be detected with a measur- 
ing system. 


Appendix 3.1: Soil Sample Preparation Prior to Gamma 


Analysis 


e Weigh the total wet soil sample in grams to 2 decimal places. Record weight as 
total wet weight. In case of a soil core, section core into desired slice increments 
and further follow the same procedure as previously described. 

e Freeze dry each soil sample (or oven dry at 105 °C). 

e Weight total dry sample (g) to 2 decimal places. Record weight as total dry weight. 


mixed and disaggregated. 


Homogenize each sample with a pestle and mortar to ensure soil material is fully 


Choose the container size according to the available sample volume to obtain a 


full (or near full) container after sample packing. Pack a suitable pre-weighed 
container (Fig. 3.2) with a subsample of the homogenised soil material: 


3 ANSI (1999). 
4ISO 11929 (2010). 
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Weigh empty container and its lid (g) to 2 decimal places; 

Pack container with soil and assure that (1) the material is firmly patted down 
inside the container, so that no void spaces are left within the sample container; 
(2) the soil surface is flat; 

Weigh packed container and its lid (g) to 2 decimal places; 

Record total packed subsample weight (g) to 2 decimal places; 

Seal container with gas-tight insulating tape; 

Label container with the pre-arranged sample ID and note date of packing; 


op 


J monao 


e Record all weighed data and sample ID onto supplied spreadsheet. 


Appendix 3.2: Preparation of a Secondary Soil Calibration 
Standard 


(A) Opening and Diluting the Mixed Traceable Standard Solution 


e Amixed traceable standard source is typically purchased from an accredited metro- 
logical institute in the form of acid liquid solution inside a flame-sealed glass 
ampoule. The vial will be opened under fume-hood conditions following the radi- 
ation protection laboratory procedures. 

Caution should be taken to avoid any spillage of the solution when the ‘swan neck’ 
is broken (along the scratch mark). 

A clean, acid-washed 100 mL Grade A glass volumetric flask is weighed (4 decimal 
places) and fitted with a small funnel. Most of the traceable reference source is 
transferred to the volumetric flask using a disposable Pasteur pipette. 

The Pasteur pipette is then filled with 1—2 mL of 4 M HCI which is used to wash the 
remaining traceable solution into the volumetric flask via the funnel. This should 
be repeated several times to ensure that all the traceable solution is transferred to 
the volumetric flask. When all the liquid has been transferred from the vial the 
liquid in the flask should be filled up to the mark with 4 M HCl and weighed again 
for 3 times. The mean arithmetic weight of the flask with the standard solution is 
recorded to 4 decimal places. 

The liquid in the volumetric flask is then transferred to a clean plastic bottle and 
labelled as ‘Stock Solution’ of the mixed radionuclide reference source and clearly 
marked with trefoil. 


(B) Dosing a Soil Sample with the Mixed Standard Solution 


e A soil sample, typically about 100 g, with the grain size <2 mm is prepared as 
mentioned in Sect. 3.3.1(i) and placed in a suitable beaker. The soil is moistened 
with a small volume of distilled water until the water level is above the soil surface, 
to enable the mixing of the slurry with the liquid standard. 

e An appropriate volume of the diluted mixed traceable reference solution is added 
to the slurry using a calibrated pipette. The amount added should be sufficient to 
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achieve <1% uncertainty of the counting statistics of the peaks of interest within 
a reasonable counting time when measuring the calibration standard. 

e The slurry is thoroughly mixed using a glass rod. 

e The slurry is dried in an oven at 105 °C or, preferably, in a freeze-drier until 
constant weight is achieved (commonly over 24—48 h). 

e The dried solid is then carefully transferred to a rotating ball mill to ensure that 
the final dried solid is completely mixed. The final product is a fine powder that 
can be allocated to appropriate sample holders, such as 50 mm diameter, dried, 
pre-weighed Petri dishes. 


For safety, any particle transfer process should be completed in a fume hood and 
the operator should be protected with a mask. 


e The filled Petri dish, and any other containers, should be re-weighed and sample 
weights recorded so the activity concentrations of each of the radionuclides in the 
standard can be evaluated. Afterwards, the standard should be bound with adhesive 
tape in distinct colour (e.g. red) for safety reasons. Its sample number, or reference 
data, should be written on the lid in indelible ink. 

e The Petri dish should then be wiped with a moist cloth to remove any extraneous 
particles. It should be kept in a sealable plastic bag when out of use. 

e These will be further used as secondary standards, in different geometries, for 
the efficiency calibration of the gamma detectors. A record should be kept on the 
radionuclide inventory: (i) present in the secondary standard, and (ii) left in the 
‘Stock Solution’. 
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Chapter 4 A) 
Conversion of Be-7 Activity get 
Concentrations into Soil and Sediment 
Redistribution Amounts 


W. H. Blake, A. Taylor, A. R. Iurian, G. E. Millward and L. Mabit 


4.1 Implementing the Event-Scale Profile Distribution 
Model 


The Profile Distribution Model (PDM) is described in detail by Walling and He (1999) 
and Blake et al. (1999). Here, alongside an overview of the basic components, we 
provide step-by-step guidance to its implementation and simple encoding within 
Microsoft Excel. 

This simple conversion model relies on exponential decline in "Be activity con- 
centration with depth wherein the shape of the profile is described by hg (See Chap. 1, 
Sect. 1.2). A first necessary step is to establish a specific ho value for the site under 
investigation. This can be achieved by plotting the depth profile from your reference 
area as an x-y plot and fitting an exponential function (see Table 4.1 and Fig. 4.1). 
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Table 4.1 Example of the depth profile data used in Fig. 4.1 where the Be concentrations have 
been determined with a 2-sigma error 


Nominal true | Total dry Mass depth Cumulative 7Be activity ae (2 
depth (mm) mass (kg m~?) mass depth concentration | sigma)? 
(<2 mm) (g) (kg m?) (Bq kg~!) 
2 83.84 0.93 0.93 70.02 4.5 
4 76.17 0.85 1.78 61.85 5.0 
6 90.14 1.00 2.78 46.91 5.2 
8 84.71 0.94 3.72 52.51 52 
10 165.91 1.84 5.56 48.06 4.4 
12 156.87 1.74 7.31 48.62 6.2 
14 175.54 1.95 9.26 24.89 5.0 
16 204.27 2.27 11.53 19.52 5.7 
18 162.85 1.81 13.34 11.03 5.4 
20 214.66 2.39 15.72 10.13 5.2 
a Analytical uncertainty 
80 5 
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Fig. 4.1 Be activity concentration against mass depth with a fitted exponential line to derive the 


value of ho 


It is important that measures of depth are represented as mass depth i.e. mass 
of soil for the core slice divided by the area of the core with units of kg m~?. This 
accounts for changes in bulk density and is a more precise measure of depth than 
true depth. họ is derived by taking the reciprocal of the exponent of the line. In the 
case of the depth profile shown in Fig. 4.1, họ is derived from: 


1 


~~ = 7.5(kg m°) 


0.133 
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Table 4.2 An example grid cell layout with example inventory data (is analytical uncertainty) 
where coding presented in Table 4.3 is required to determine h, Cg and h’ (see Chap. 1, Appendix 
1.1), without particle size correction. Calculated values following Table 4.3 are shown in italics 


A B | c D E F | G 
Sample "Be t Uncorrected | Eroded soil Deposited Uncorrected 
ID inventory erosion "Be activity soil Be deposition 
amount (h) conc. activity conc, amount (h') 


To run the PDM in Excel, the key components of the model can be coded into a 
cell grid as illustrated in Table 4.2. In this example, the reference inventory value is 
set at 524 + 46 Bq m~? (n = 15 where uncertainty represents 2 standard deviations). 
The sample points 1—10 are erosional (i.e. an inventory deficit as compared to the 
activity recorded in the reference area) and sample points 11 and 12 are depositional 
(i.e. an inventory excess as compared to the activity recorded in the reference area). 

Coding required for columns D, E, F and G to calculate h, Cy and h' respectively 
(see Chap. 1, Appendix 1.1 for relevant model equations) is provided in Table 4.3. 

The PDM may also be applied using a specific Microsoft Excel Add- 
in—developed by Walling and He (1999) and refined by Walling et al. (2002)—which 
is available at the IAEA website (see http://www-naweb.iaea.org/nafa/swmn/models- 
tool-kits.html). 


4.2 Accounting for Size Selectivity of Erosion 
and Deposition Processes in the Profile Distribution 
Model 


As pointed out in the Sect. 2.4 of the Chap. 2, some advanced users will need to 
account for particle size selectivity associated with soil erosion processes to convert 
inventory loss into soil erosion amounts. Different approaches are described in the 
literature (see Taylor et al. 2014 or Yang et al. 2013). Here, we describe the method 
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Table 4.3 Coding for columns D, E, F and G in Table 4.2 to derive the key elements of the PDM 
as described in Sect. 1.2 


Cell number Code Explanation 

(assuming 

layout as per 

Table 4.2) 

D2 =$I$2 x LN($I$3/B2) (see This Excel code is to derive ‘h’ i.e. the 
Chap. 1 Eq. 1.3) mass depth of soil lost (NB as positive 


value) i.e. Aref/A to account for the 
observed inventory deficit in cell D2 
where cell I2 contains the value for ho 
(in this example 7.5 kg m=”) and cell 13 
contains the reference inventory (in 
this case Aef = 524 Bq m=). Note that 
this code can then be copied down to all 
eroding site points 


E2 =($1$3 — B2)/D2 (see This Excel code permits to estimate the 
Chap. 1 Eq. 1.5) 7Be concentration of the soil eroded from 
this sample point or zone 
F11 =(D2 x E2)+...+(D11 x E11); | This Excel code allows to calculate the 
SUM(D2:D11) (see mean activity concentration of upslope 
Chap. 1 Eq. 1.6) eroded soil used in the equation for h’ 


below. Note that all upslope contributing 
points need to be included in the code 
across (shown in shortened form). This 
needs to be copied for all depositional 


points 
G12 =(B12 — $I$3)/F12 (see This is to calculate h’ i.e. the mass depth 
Chap. 1 Eq. 1.4) of soil gained to account for the inventory 


excess in cell D11. This needs to be 
copied for all depositional points 


proposed by Taylor et al. (2014) and compare to other methods so users can make 
an informed decision about which approach to take. 

The PDM assumes that the soil profile is top-sliced by erosion and that loss of 
inventory is directly related to the mass depth of soil lost (see Fig. 1.3). However, 
enrichment of ‘Be in fine soil fractions (Fig. 4.2) means that these can be prefer- 
entially removed by erosion processes resulting in an augmented inventory deficit 
relative to erosion amount (Blake et al. 2009; Taylor et al. 2014). It should be noted 
that in some cases where soil has an overall fine texture and rainfall is of high mag- 
nitude, particle size correction may be less important (e.g. Porto and Walling 2014) 

The approach suggested by Taylor et al. (2014) accounting for particle size selec- 
tivity follows the principles of the method proposed by He and Walling (1996) 
wherein particle size correction factors P and P’ for eroding and depositional sites 
respectively are calculated based on measurements of ‘geometric’ Specific Surface 
Area (SSA unit in m? g7!). P is derived from: 
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Fig. 4.2 Enrichment of 805 
stable Be, as a surrogate for 70 | 
7Be, in two fine grained soil 
materials with contrasting 
texture where soil A has a 
greater sand content (Taylor 
et al. 2014) 


y = 52.21 1x071% 
R? = 0.99 


a Soil A 
© Soil B 


Be (mg kg') 
A 
o 
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y = 16.623x 94343 
R? =0.97 


r r T T j 
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Specific Surface Area (m? g`!) 


Sm E 
P = | — 4.1 
( 3 ) a 
where Sm is the mean SSA of the mobilised fraction and S, is the mean SSA of 


the bulk soil material and v is a constant describing the relationship (i.e. the power 
function exponent) between 7Be concentration and SSA (Fig. 4.2). 


P’ is derived from: 
Sa\° 
P' = | — 4.2 
( Sm ) ( ) 


where S4 is the mean SSA of the deposited material. 

The SSA of soil and sediment samples is most commonly performed using laser 
granulometry. It is also possible to derive SSA from particle size data derived by grav- 
itational sedimentation methods, most commonly the pipette method (Day 1965). 

It is important at this stage to consider that amounts of erosion estimated using 
a conversion model consider inventory deficit i.e. the amount of Be inventory lost 
from a specific area due to soil erosion. This inventory deficit is converted to mass 
depth soil loss using the exponential shape of the ‘Be distribution in the soil profile 
(Eq. 1.1). In this context, Taylor et al. (2014) propose that when correcting erosion 
rates (mass depth lost, kg m~?) for particle size enrichment it is, therefore, necessary 
to apply the correction factor to the inventory deficit only (i.e. A minus A;er) which 
has the effect of suppressing the loss in inventory according to the particle size of the 
fraction mobilised. If fine sediment is preferentially mobilised by erosion processes, 
failure to account for this will lead to overestimation of erosion. 

In the case of the above, following the principles of He et al. (2002), a particle-size 
corrected erosion rate, h,, can be calculated using the following equation (compare 
to Eq. 1.2): 


Aref 


he = hol 
= Aref = (Aref aa A)/P 


(4.3) 
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For estimates of deposition, the 7Be concentration of deposited material needs 
to be adjusted according to the particle size of the deposited fraction relative to the 
eroded fraction. The equation to provide corrected deposition rates, h’, is (compare 
to Eq. 1.4): 


where C4c is the weighted mean concentration of sediment mobilised from the ups- 
lope contributing area (Blake et al. 1999; Walling et al. 1999) which has been cor- 
rected for particle size enrichment during erosion, and then corrected for deple- 
tion during deposition (He et al. 2002). The particle size corrected concentration of 
deposited material from an individual eroding point (ca’) can be estimated using the 
following function: 


cd’ = P P'A,ep(1 = el) 1h (4.5) 


These data processing requirements need to be considered in the context of the 
sampling challenges detailed in Chap. 2. 

The above equations can be used to modify the Excel-based PDM model coding 
provided in Table 4.3. 

There are several challenges with particle size correction that all users must be 
aware of. Determining the relationship between the activity concentration and SSA 
to gain the value of the exponent, v, in Eqs. (4.1) and (4.2), requires particle size 
separation experiments (Taylor et al. 2014) and use of a laser granulometer or quan- 
titative sieving (and application of SSA theory which makes assumptions about 
particle sphericity). Sampling for S, and Sg is relatively straightforward (sampled in 
field from stable and deposition sites) but sampling for Sm is less straightforward. 
Options for this include rainfall simulation experiments on the hillslope after the ero- 
sion event to collect representative eroded material, installation of Gerlach troughs in 
study field prior to anticipated erosion events or estimating SSA of eroded material 
from remaining soil-and comparison with uneroded material (see Yang et al. 2013) 

All particle size corrections require a necessary simplification of reality since soil 
characteristics are highly variable within and between sites. Example data illustrating 
the influence of particle size correction approaches on soil erosion estimate results 
are given in Sect. 4.4. 


4.3 Extended Time Series Conversion Model: The Theory 


If an advanced user wishes to apply the ’Be approach over a longer time period, for 
example over a wet season that contains notable erosional events in series, an alterna- 
tive approach is required (Walling et al. 2009). Modifications to the PDM approach 
are required to account for inventory flux (Bq m7”), decay and the erosivity of events. 
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Fig. 4.3 Relationship 605 _— 
between the percentage F] ee a ae 
under-estimation of the 50 | ao 

erosion rate versus the length 40 a 


of the study period (from 
Walling et al. 2009) 
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Failure to account for these factors has been shown to lead to underestimation of 
redistribution rates (Walling et al. 2009). In this approach, the cumulative change 
in inventory at the sample site is considered rather than a direct comparison to the 
stable reference value for each day. A reference inventory for the field site is required, 
however, to provide an indication of erosion or deposition at each individual sample 
location. 

The basic principle of the extended time-series model is that the amount of erosion 
estimated by ’Be inventory deficit and deposition estimated by ’Be inventory excess 
needs to be apportioned across the time period of study while allowing for the effects 
of radioactive decay. The problem of the PDM over extended time periods is con- 
sidered by Walling et al. (2009) for a hypothetical scenario where: (1) the reference 
inventory begins at zero due to extended dry period, (2) a single rainfall event of same 
intensity takes place each week with same erosivity, (3) ho remains the same. The 
results (Fig. 4.3) demonstrate that the 7Be approach increasingly underestimates the 
erosion rate with substantial error as time progresses into the study period (e.g. 30% 
after 2-3 weeks; 50% after 3 months). In addition to this issue of decay of the ear- 
lier erosion event signals, there are other potential complications linked to temporal 
distribution of ’Be during the study period, temporal distribution of erosion (rainfall 
erosivity), potential variability of other key parameters through time, especially on 
freshly ploughed surface, e.g. ho. 

In this context, every user needs to consider the sensitivity of their study data to 
these issues and either take action or acknowledge these limitations and the uncer- 
tainty introduced into their results. 

In this section we present a step-by-step approach to implementing the extended 
time series model in Microsoft Excel based on Walling et al. (2009). There are three 
main components: (1) calculating relative erosivity of rainfall, (2) estimating soil 
loss and (3) estimating deposition. 

The aim of the extended time series model is to allow 7Be measurements to be 
extended over a longer timeframe that is more representative for some environments 
for example where erosion takes place over a prolonged wet season. 

In addition to the evidence collected for the PDM, the extended time series model 
requires a reconstructed record of daily fallout deposition (Bq m~? d7!) for the 
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whole study period and a time series of daily (where appropriate) relative erosivity 
(see below). Together, these account for the effects of time variant erosion and the 
time variant fallout as well as for the radioactive decay. 

The first step in applying the model is the same as for the PDM (Sect. 4.1) namely 
to determine whether the site is erosional or depositional i.e. compare inventory at 
erosional point (A) or depositional point (A’) to reference inventory Aef- 

The second step is to build a mass balance model framework used to define the 
7Be inventory at the end of each day A(t) (Bq m7). This is determined from: (i) the 
7Be inventory of the previous day [A(t — 1)], (ii) effects of radioactive decay [A], 
(iii) any fallout input (that day [F(t)], (iv) any loss in "Be inventory due to erosion 
[Aioss]. Of course, the mass balance must ultimately be solved for the latter element. 
These components are linked together as follows: 


A(t) = A(t — 1) exp(—A) + F(t) — Aloss (1) (4.6) 


Ajoss (t) in Eq. (4.6) will reflect (4) amount of soil eroded [R(t)] as a mass depth 
i.e. kg m7”, (ii) existing inventory at the point [A(t)], and (iii) depth distribution of 
inventory i.e. họ (for that day) i.e.: 


R(t) 
Aloss(t) = [A(t — 1) exp(—A) + F (t)] * [i = exp(- 7) (4.7) 


In the prior equations, the amount of soil eroded R(t) is assumed to be proportional 
to the rainfall erosivity for that daily time-step and can be expressed as following: 


R(t) = E,(t)* C (4.8) 


In Eq. (4.8), C is a constant that links relative erosivity (i.e. E,) to erosion amount. 
The value of C differs between sample points but at any one sample point it remains 
constant through the study period. 

To derive C and hence R(t), the user must establish a continuous mass balance for 
the study period for each sampling point experiencing erosion. Guidance on setting 
this up in Microsoft Excel is given in the next section. The value of A(t) at the end of 
the study can then be related to the measured inventory (A) at the time of sampling. 
Note the value at the start is equal to the measured reference inventory at the beginning 
of the study (often zero due to cultivation of extended dry period). Given all the data 
that have been derived up to this point, the mass balance can be solved for C using 
the optimising Solver routine in Microsoft Excel. This so called ‘what-if analysis’ 
tool basically finds the optimum value of C to explain the measured inventory within 
the constraints of the mass balance detailed in the worksheet. 

Application of C to the E,(t) values in the time series then gives the estimated 
erosion amount for the sample point. An example of spreadsheet coding is given in 
Sect. 4.4. 
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As noted above, in the extended time-series model, a key part of the process is 
using rainfall erosivity to apportion erosion potential across the study period. For the 
study period in question, a measure of relative erosivity is required for each event 
occurring during the study period. 

Arguably, we need knowledge of site specific thresholds for erosion initiation 
based on rainfall patterns and dynamics and site-specific conditions need considera- 
tion e.g. runoff initiation processes. Walling et al. (2009) used high resolution rainfall 
data (minimum 30 min interval data) to evaluate temporal distribution of erosivity 
(E,). This approach is linked to RUSLE daily values of erosivity based on the product 
of the total kinetic energy of rainfall (E) and maximum 30 min rainfall intensity (739) 
i.e. E39. This means a continuous rainfall record at study site is required (30 min 
totals). It should be noted that a measure of relative erosivity in this context is more 
important than an absolute value. Local empirical data, if available, may be more 
appropriate at some sites. Taylor (2012) followed procedures outlined by Yin et al. 
(2007) and Morgan (2001) to derive EIo following Schuller et al. (2010). 


4.4 Implementing the Extended Time Series Model 


The extended time series model can be implemented in a variety of ways. We propose 
here one version for coding in Microsoft Excel as an example to develop from. Such an 
approach permits the user to see the model working step-by-step and to troubleshoot 
any coding errors or data anomalies that emerge. 

The model structure comprises a series of columns that represent the time series 
of daily fallout, erosion etc. (Table 4.4) and a set of ‘model interface’ cells where 
primary data is inputted and results are calculated. 

To permit cell O7 to show the most likely erosion rate, the Solver Routine must be 
run. This requires the actual measured inventory at the sample point to be entered into 
the Solver dialogue box (Fig. 4.4) and the routine run so the value of O5 is adjusted 
to equal the measured inventory. The routine does this by adjusting the unknown 
constant C until the value of O7 is equal to the inventory measured. 

The extended time series model is an appropriate approach to address the challenge 
of using ’Be across multiple rainfall events. In addition to the requirements of the 
PDM, it requires rainfall data and site-specific determination of relative erosivity. 
Implementation of this advanced model requires careful coding and usage of data. 
Users need to be aware of specific sample point attributes and relationships to other 
points and to ensure that each working step is traceable in case of a need to revisit data. 
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Table 4.4 Microsoft Excel coding for time series and mass balance columns. It is assumed user 
has a basic knowledge of Microsoft Excel and that the column headings are entered in the first row 


Cell number Code Column title and explanation 

A2 (and relevant | n/a Date of time step (in DD/MM/YYYY format). 

cells below) This should be copied down for the duration of 
the study period from the start of rainfall 
measurement to the final erosion event to be 
studied 

B2 (and relevant | n/a Daily inventory fallout flux (F) value 

cells below) (Bq m7) which has been derived from a 
rainfall measurement and sampling programme 
described in Chap. 2 

C2 (and relevant | =0.0130 Decay constant for 7Be for 1 day time step 

cells below) 

D2 n/a The study area initial plot reference inventory 
at the start of the rainfall measurement and 
monitoring, typically zero (Bq m~?) 

D3 =B3 + (D2 * EXP(—C3 x 1)) Daily reference inventory for time step based 


on decay corrected inventory from prior day 
plus the fallout input of the current day 

(Bq m™?). Note this column is copied down to 
the end of the study period 


E2 (and relevant | n/a Relative Erosivity (E+) value for the timestep. 

cells below) Note that all E, values in the column will sum 
to 1. Note this column is copied down to the 
end of the study period 

F2 =$0$2 Constant C. Note this code links it to a cell in 
the model ‘interface’ cells. Note this column is 
copied down to the end of the study period 

G2 =E2 x F2 Relative Erosion Rate (R). Note this column is 
copied down to the end of the study period 

H2 (and relevant | =G2/$0$3 Particle size corrected relative erosion rate 


cells below) 


where cell O3 contains the simple correction 
factor P in the model ‘interface’ cells area 
(based on Yang et al. 2013). See Sect. 4.2. Note 
this column is copied down to the end of the 
study period 


I2 


=(D2) x (1 — EXP(G2/$0$4)) 


Inventory loss (Ajoss) (Bq m2) where cell 04 
contains the value of hg. This is for first time 
step. For following time steps, cell I3 then 
needs updating as below to account for decay 
or previous time-step 


13 


=(J2 x EXP(—C3) + B3) * (1 — 
EXP(G3/$0$4)) 


Note this column is copied down to the end of 
the study period 


J2 


=D2 — I2 


Modelled inventory (A) (Bq m7“) for first time 
step. For following time steps, see cell J3 as 
below to account for decay or previous 
time-step 


J3 


=J2 * EXP(—C3 x 1) + B3 — I3 


Note this column is copied down to the end of 
the study period 


K2 


=IF(G2 <>0, 12/G2, 0) x $0$3 


Activity of eroded soil (Ce) (Bq kg7!). This is 
used in deposition calculations downslope 
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E 
Solver Parameters =S 


| Set Objective: soss =) 
To: Max Min © value OF: ENTER MEASURED INVENTORY 


By Changing Variable Cells: 
$Os2 


E 


Subject to the Constraints: 


Add 

Change 
Delete 

| Eo 
Reset All 

| 

| j| Load/Save 

Make Unconstrained Variables Non-Negative 
Select a Solving GRG Nonlinear [z] Options 
| Method: 


Solving Method 


Select the GRG Nonlinear engine for Solver Problems that are smooth nonlinear. Select the LP 
Simplex engine for linear Solver Problems, and select the Evolutionary engine for Solver 
problems that are non-smooth. 


| | Help | Solve Close 


£ d 


Fig. 4.4 Dialogue box for Miscrosoft Excel Solver command set up to work with the user interface 
cells in Table 4.5. Note where the value of the measured inventory needs to be entered 


4.5 Hillslope Sediment Budget Examples and Inclusion 
of Uncertainty 


The data presented in Table 4.2 can be used to check model coding by comparing 
model outputs derived with this dataset to the results presented in Table 4.6. It is 
important to note that Table 4.6 is provided for checking of the central result only 
and does not include an uncertainty. It is critical that users report their soil erosion 
estimates with an appropriate uncertainty and that the means by which the uncer- 
tainty has been determined is transparent to the reader of reports and subsequent 
publications. 

Uncertainty can be derived from a range of different factors within the sampling 
and analysis process (Fig. 4.5). The most obvious source of uncertainty is the ana- 
lytical error originating from the gamma spectrometry measurement (see Chap. 3) 
which can easily be propagated through to the resulting soil erosion estimate. This 
is, however, arguably less environmentally relevant than other factors. An alternative 
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Table 4.5 Microsoft Excel coding for user interface cells (linked to cells described in Table 4.4). 
Note that until the Microsoft Excel Solver is run using the coding provided below, the numbers in 
the results cell are arbitrary. Using the ‘Solver’ command the modelled inventory can be fitted to the 
measured inventory to allow cell OS to equal the sample point inventory by changing the constant. 
Values of ho and particle size correction can be changed manually in cells 03-4 for scenario testing 


Cell Code Column title and explanation 
number 
02 Enter any start value e.g. 1 Cell where value of C will be calculated 
by the Solver routine 
03 n/a Enter particle size correction factor P 
(enter 1 if no correction required) 
04 n/a Enter value of ho 
05 Enter the cell code for column J at the This cell value will change depending on 
last date in the time series e.g. = J128 the value of the constant C. Initially the 
would be appropriate for a record witha | number showing will be an artefact of 
127 day extended time series period the entry in cell O2 until the solver 
routine has been run 
07 =SUM(H2:H##) where ## is the row Total erosion based on sum of all values 


number corresponding to the last time 
step in the time series 


of R (kg m2). This is the final result but 
should only be recorded after the solver 
routing has been run for the worksheet 
for the specific sample point 


Table 4.6 Model output data using the demonstration dataset in Table 4.2 (where P = 1.61 and P’ 


= 0.82) 

Sample ID PDM soil redistribution | PDM soil redistribution | PDM soil redistribution 
amount (no particle size | amount (particle size amount (particle size 
correction) (kg m~*) corrected Taylor et al. corrected Yang et al. 

2014) (kg m~?) 2013) (kg m~?) 

1 =79 —3.9 —4.9 

2 =2.1 -1.3 =13 

3 —4.2 =23 —2.6 

4 —6.7 —3.4 —4.1 

5 —4.1 =23 —2.5 

6 —4.7 —2.6 —2.9 

I —2.8 —1.6 =1.7 

8 =2.1 —1.2 =1.3 

9 —3.5 —2.0 =2.2 

10 —1.4 —0.8 —0.9 

11 5.7 4.5 3.2 

12 10.5 8.4 6.0 
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in fallout 
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Spatial variability in 
within reference PS depth profile (ho) 


Spatial variability within i Gamma counting 
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Particle size and 
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Preferential mobility 
(physical and geochemical) 


Fig. 4.5 Potential sources of uncertainty in 7Be-based estimates of soil erosion (Taylor et al. 2013) 


and more relevant approach is to work with the uncertainty derived from repeated 
reference inventory measurements. Herein, the coefficient of variation between cores 
represents a real-world measure of local variability which should be accounted for 
in erosion estimates. In practice, the model can be run with the upper and lower 
limits of the reference inventory and the range of outputs used accordingly for each 
sample point. It should be noted that results from points that differ marginally from 
the inventory might switch between eroding and depositing depending in reference 
uncertainty. Similarly, where a spatial integration approach (see Chap. 2) has been 
adopted in plot sampling, variability in inventory within zones can be considered. 
Whichever approach is selected, the user must report what was performed with full 
transparency. 

Any approach taken to determine and report uncertainty will have inevitable site- 
specific and study-specific caveats. Users are encouraged to explore potential sources 
of uncertainty and undertake suitable sensitivity analysis when running conversion 
models to convince the scientific and stakeholder communities of data quality and 
the environmental relevance of results. 


Glossary 


Profile Distribution Model (PDM) Simple model to convert 7Be inventory loss 
into soil erosion amounts by linking inventory 
change to the depth profile to derive the mass 
depth of soil loss. This model is applied to sin- 
gle events of short period of time. 

Extended time series model Improved conversion model that permits the 
effect of radioactive decay on soil erosion rate 
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estimates using Be over a longer period (e.g. 
wet season) to be accounted for. 

Sensitivity analysis Process to link uncertainty in the output of a 
model to different sources of uncertainty in its 
inputs with the aim of an increased understand- 
ing of the relationships between input and output 
variables. 
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Chapter 5 A) 
Research into Practice—Linking Be-7 get 
Evidence to Land Management Policy 

Change for Improved Food Security 


M. Benmansour, W. H. Blake and L. Mabit 


5.1 The Importance of the Short-Term Perspective 
for Land Management Policy Makers 


The short-term perspective on soil redistribution in the landscape offered by Be 
is essential to policy makers with responsibility for developing efficient land man- 
agement strategies to support food and water security. Soil conservation is vital for 
enhancement of food production on hillslopes of agroecosystems. Mobilisation of 
eroded soil and transfer downstream leads to siltation of river channels, lakes and 
reservoirs, which presents a credible threat to river basin ecosystem service provi- 
sion and water security. In addition, energy security is threatened by siltation behind 
hydropower dams. The short-term perspective of this cosmogenic radioisotopic tool 
means it can provide relatively rapid assessment of very recent changes in practice. 
Because of its agro-environmental behaviour, this short-lived radioisotope a reliable 
natural isotopic tracer to assess the effectiveness of recent soil conservation strategies 
(Mabit et al. 2008; Taylor et al. 2013). 

The key challenge that many policy organisations face is obtaining convincing 
and scientifically sound evidence upon which realistic and wise informed decisions 
can be based. While many conventional methods are capable of determining on site 
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erosion rates at a point (e.g. rainfall simulation, erosion pin measurements, survey 
techniques), the scaling of these measurements to the wider hillslope or landscape 
unit is challenging. This is especially the case given that equipment associated with 
many conventional approaches interferes with either farming activities themselves 
or natural soil erosion processes. 

An advantage of "Be is that this isotopic tracer is integrated into the natural hydro- 
logical cycle. Delivered naturally via rainfall to the soil surface and redistributed via 
erosion processes, it operates without any interference by deployment of measure- 
ment equipment. In addition, the method can be applied retrospectively and sampling 
strategies driven by field observations after an erosion event. The method is not with- 
out complications and challenges but if undertaken with care and due attention to 
assumptions and limitations (as described in Chaps. 1—4), the above advantages offer 
opportunity to gain critical information on soil and sediment redistribution to support 
policy decisions on soil conservation and environmental management. 


5.2 Linking Nuclear Techniques in Soil Erosion 
and Conservation to Policy Change 


Early studies conducted in the 1990s outlined the principles and opportunity for use 
of Be as a soil erosion and soil redistribution tracer (Wallbrink and Murray 1996; 
Blake et al. 1999; Walling et al. 1999; Blake et al. 2002). This foundation-laying 
research has since had a wider impact on soil conservation and environmental man- 
agement policy through applied studies of the method around the world. This has been 
facilitated by research development and refinement performed within [AEA’s Coor- 
dinated Research Projects (CRPs) such as the CRP D1.50.08 “Assess the effectiveness 
of soil conservation techniques for sustainable watershed management using fallout 
radionuclides” (2002-2008), the CRP D1.20.11 “Integrated Isotopic Approaches 
for an Area-wide Precision Conservation to Control the Impacts of Agricultural 
Practices on Land Degradation and Soil Erosion” (2009-2013) and more recently 
through the on-going CRP D1.50.17 “Nuclear Techniques for a Better Understand- 
ing of the Impact of Climate Change on Soil Erosion in Upland Agro-ecosystems” 
(2016-2021). Following the CRPs activities, several guidelines were produced by the 
SWMCN Subprogramme of the Joint FAO/IAEA Division of Nuclear Techniques in 
Food and Agriculture (e.g. IAEA 2014; Mabit et al. 2018). 

Impact on policy has been further extended worldwide to several IAEA Member 
States and FAO Member Countries by methodology transfer and capacity-building 
activities within [AEA national and/or regional technical cooperation projects (e.g. 
RAF5075 “Enhancing Regional Capacities for Assessing Soil Erosion and the Effi- 
ciency of Agricultural Soil Conservation Strategies through Fallout Radionuclides 
in Africa”) which includes targeted training courses on the use of Be as a tracer 
to support soil conservation policy by leading researchers. As a result of these out- 
reach and science-policy impact activities, application of this specific methodology 
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has been extended to a range of soil conservation challenges, in collaboration with 
national and international organisations responsible for environmental management 
policy to support food and water security or management of forest resources. 

For example, Blake et al. (2009) undertook a landscape sediment budget approach 
with "Be alongside !*’Cs and *!°Pb to evaluate the post-fire loss of topsoil and nutri- 
ents from forest affected by wildfire to a major water supply reservoir in Australia, 
in collaboration with the government agency responsible for developing policy on 
fire management for water quality and quantity. Under future climates, wildfire will 
exacerbate threats to water security. The applied fallout radionuclide (FRN) research 
demonstrated that burning of surface vegetation can invigorate hillslope hydrological 
response with marked increases in sediment and nutrient delivery to river networks 
and reservoirs. Negative water quality effects include for example high turbidity, 
toxic algal blooms and fish kills with implications for water supply at critical times 
in the water year. Through quantifying post-fire runoff and nutrient yield processes, 
the research informed catchment management decisions, policies and water resource 
risk assessment in Australia. Similarly, researchers in Chile, applied ’Be using a tran- 
sect approach (Schuller et al. 2006) for evaluation of the impact of forest harvesting 
on soil erosion fine sediment yields and for contribution to improve the effective- 
ness of mitigation measures before, and after clear-cut operations and during the 
establishment of the new plantations. This work demonstrated the benefit of using 
7Be to document the effects of recent forest clearance upon slope soils in a timber 
harvesting operation. The short-term perspective of "Be meant the use of mitigation 
strategies could be assessed and, in this case, a high sediment delivery ratio indicated 
to policymakers the need for further soil conservation strategies. 


5.3 Example Impact Case Study: Support Provided by The 
7Be Technique to Shape Soil Conservation Policy 
in Morocco 


5.3.1 Soil Erosion and Conservation Policy Challenges 
in Morocco 


Soil erosion is the main land degradation process in Morocco, which is affecting 
up to 40% of its land area. About two million hectares of Moroccan agricultural 
lands are affected by soil erosion by water (Dahan et al. 2012). On average, soil 
erosion ranges from 5 to 20 t ha~! yr~!, but exceeds these magnitudes in north- 
ern and north-western basins. For example, in the pre-Rif hills region, soil erosion 
exceeds 50 t ha! yr~'(HCEFLCD 2017). Out of a total area of 20 million ha of 
watershed in Morocco, 50% is estimated to be subjected to very high erosion risks 
with a yearly soil loss of around 100 million tons. This annual soil loss leads to a 
reduction of 75 million mê of downstream dam water storage capacity. Each year 
about 0.5% of the country’s reservoir storage capacity is lost. These soil erosion 
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and conservation challenges require both spatially representative data on soil ero- 
sion amounts and sediment delivery ratios. The highly dynamic nature of land use 
change and implementation of new conservation policies, however, means that long- 
or medium-term FRN approaches are not recommended to deliver the data required 
by policymakers on the extent of the problem and the level of success in conservation 
approaches applied (Mabit et al. 2008). Consequently, researchers in Morocco have 
been using the ’Be methodologies outlined in this book to tackle these challenges 
and the obtained data have been used to assist development of new soil conservation 
policy in a wide range of settings. 


5.3.2 From the Laboratory to the Field: The Approach 


FRN methodologies have been successfully used by the “Centre National de 
l Energie, des Sciences et des Techniques Nucléaires” (CNESTEN) in many cul- 
tivated areas of Morocco for supporting soil conservation strategies tested or imple- 
mented by the government. The research has shaped soil conservation policy with 
the collaboration of “Institut National de Recherche Agronomique (INRA) from the 
Ministry of Agriculture and Maritime Fisheries and “Centre de Recherche Forestière” 
(CRF) from the “Haut Commissariat des Eaux et la Lutte contre la Désertification” 
(HCEFLCD). Fallout "Be was applied in combination with '°’Cs and !°Pb to esti- 
mate soil erosion in agricultural or forest fields under different land uses in the regions 
of Rabat, Tétouan, Casablanca and Fes (Fig. 5.1) both before and after implementa- 
tion of recommended soil conservation approaches. The transect sampling strategy 
was generally adopted (see Chap. 2) and FRNs were determined at CNESTEN Lab- 
oratory using gamma detectors (see Chap. 3). Key studies where the "Be approach 
has been integrated into policy-development research: 


(1) In semi-arid and Mediterranean climatic areas (Rabat and Tétouan regions), a 
no-till technique was tested as a soil conservation practice and compared to 
conventional tillage in the experimental sites “Marchouch” and “ Herchane” of 
the region of Rabat and Tétouan respectively; 

(2) In semi-arid climatic areas (Casablanca-Settat region), Be was used to evaluate 
the efficiency of conservation practices involving Atriplex plantations (imple- 
mented based on results from other FRNs) and, cereal and fruit plantations 
effectiveness were assessed within the “Oued Mellah” watershed management 
programme; 

(3) In variable climatic areas from cold-humid to semi-arid climate (Fes region), 
the efficiency of soil conservation practices based on the use of Allepo pine trees 
combined with dry stones and gabion and also fruit plantations was assessed 
within the framework of the “Allal Fassi” watershed management programme. 
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Fig. 5.1 Locations of the study areas = (1) Marchouch site (region of Rabat), (2) Herchane site 
(region of Tétouan), (3) Oued Mellah watershed (region of Casablanca) and (4) Allal Fassi (region 
of Fes) 


5.3.3 Key Findings and Policy Impacts from ’Be Application 
in Morocco 


Previous investigations using other FRNs (!37Cs or 7!°Pb.x) allowed establishment 
of the long-term soil erosion rates of the three investigated regions, which ranged 
from 8 to 58 t ha~! yr! (Benmansour et al. 2013; Benmansour et al. 2016; Yassin 
et al. 2017). The net soil erosion rates appeared to be closely related to the rainfall, 
slope and the past land use. The sediment delivery ratio in all areas is generally high 
reaching 80 to 100% which led to implementation of soil conservation strategies to 
keep fertile soil material on the hillslopes. The short-term perspective of Be enabled 
the effectiveness of these soil conservation strategies to be assessed and underpinned 
their wider adoption in the landscape. Figure 5.2 reports the mean values of short- 
term erosion amounts estimated at different study sites and corresponding to fields 
under conventional tillage or without management and those corresponding to fields 
under soil conservation practices. 
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Soil redistribution data from 7Be measurements indicated that soil loss had been 
reduced significantly under no-till as compared to conventional tillage in Rabat and 
Tétouan regions. Indeed, soil erosion rates were lowered by 50% for the Marchouch 
site (Fig. 5.2a) and by 40% for the Harchane site (Fig. 5.2b). At the Oued Mellah 
watershed, the results highlighted that high density Atriplex plantations reduced soil 
loss by approximately 57 to 80% compared to Atriplex plantations with low density 
while for the site under fruit plantations and cereals, soil erosion has been decreased 
by 58% compared to bare soils (Fig. 5.2c). For Allal Fassi Watershed, erosion was 
lowered by 54% for soil under young Allepo pine, dry stones and gabions compared 
to soil under only Allepo pine and dry stones and by 51% for agricultural fields under 
cereals and fruit plantation as compared to bare soils (Fig. 5.2d). These findings have 
emphasised the need for wider adoption of these proven soil conservation approaches 
in the Moroccan landscape. Indeed, the no-till practice tested by Moroccan INRA 
researchers in their experimental sites is being applied and transferred by the Ministry 
of Agriculture to agricultural farms that have similar agro-environmental conditions 
as the experimental stations of INRA. The demonstrated high performance of soil 
conservation practices tested by HCEFLCD in the investigated watersheds means 
they will be adopted at a wide scale. 

The relevant results derived from the application of these FRNs for short and 
longer-term erosion rates in Morocco were presented jointly by CNESTEN and 
HCEFLCD at the 12th Session of United Nation Convention for Combating Deser- 
tification, COP 12, Turkey, October 2015. The impact of this research on policy was 
emphasised by a lead policy maker, the Head of the Forestry Research Centre of 
Morocco’s High Commission of Forest and Water and Combating Desertification 
(HCEFLCD). He spoke about the role of isotopic techniques in tackling the effects 
of climate change, particularly drought leading to a reduction of up to 75% in grain 
yields: “Using isotopic techniques, we were able to accurately assess soil erosion 
and the effectiveness of soil conservation practices and make concrete recommen- 
dations to policy makers... These all delivered real change for people who rely on 
the land for their livelihoods.” (IAEA 2016). 


5.4 Future Trends and Opportunities in Using Be 


The above science-to-policy examples in Morocco illustrate the fundamental role that 
the short-term perspective on soil redistribution rates provided by "Be can play in 
the development and implementation of soil erosion policy. The need for such infor- 
mation has never been so crucial. Against threats from soil erosion, it is predicted 
that global food production must increase by 70% to feed the projected growth of 
the world’s population (from 7 to 9 billion by 2050) (FAO 2017). Extreme weather 
(wet and/or dry) events may become more frequent with changing global climate 
with impacts on soil quality and erodibility. Developing community resilience to 
such recurrent climatic events demands “win—win’ systems of soil conservation that 
permit enhancement of land productivity while protecting ecosystem services in the 
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(a) Machouch Site (Rabat region) (b) Herchane site (Tétouan region) 
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Fig. 5.2 Short soil erosion rates under different land uses associated with the study sites: Marchouch 


(a), Herchane (b), OuedMellah (c) and AllalFassi (d). LD Low density; HD High density; E East 
exposure; W West exposure 
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wider river basin. Key information and guidance provided by nuclear and isotopic 
tools such as the ’Be technique provides valuable support to convince policy mak- 
ers and farmers to adopt and promote widely existing effective climate smart soil 
conservation strategies. 
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